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Abstract. Some types of bacteria use rotating helical flagella to swim. The motion of such organisms takes 
place in the regime of low Reynolds numbers where viscous effects dominate and where the dynamics is 
governed by hydrodynamic interactions. Typically, rotating flagella form bundles, which means that their 
rotation is synchronized. The aim of this study is to investigate whether hydrodynamic interactions can 
be at the origin of such a bundling and synchronization. We consider two stiff helices that are modelled 
by rigidly connected beads, neglecting any elastic deformations. They are driven by constant and equal 
torques, and they are fixed in space by anchoring their terminal beads in harmonic traps. We observe that, 
for finite trap strength, hydrodynamic interactions do indeed synchronize the helix rotations. The speed 
of phase synchronization decreases with increasing trap stiffness. In the limit of infinite trap stiffness, the 
speed is zero and the helices do not synchronize. 

PACS. 05.45.Xt Synchronization; coupled oscillations - 47.15.Gf Low-Reynolds-number (creeping) flows 
- 87.16.Qp Pseudopods, lamellipods, cilia, and flagella - 87. 19. St Movement and locomotion 



1 Introduction 

Many types of bacteria, such as certain strains of Es- 
cherichia coli or Salmonella typhimurium, swim by rotat- 
ing flagellar filaments, which are several micrometers long 
and about 20 nm in diameter (the size of the cell body 
is about 1 fim) [1-5]. The complete flagellum consists of 
three parts: the basal body which is a reversible rotary 
motor embedded in the cell wall, the helical filament that 
acts as propellor, and in-between a short flexible coupling 
called the proximal hook [2-5]. The motor is powered by 
protons moving down an electrochemical gradient [2,4,5], 
which generates a constant torque independent of the dy- 
namic load [4,6]. The rotation rates for the flagella of 
freely moving bacteria are of the order of 100 Hz [3-5]. 
The filaments are polymers with high fiexural and tor- 
sional stiffness [2,4,5]. However, they are flexible enough 
to switch between different helical forms with distinct cur- 
vature and twist [3-5]. 

Typically, the filaments rotate in synchrony, i.e., the 
helices are locked in phase so that they can form bundles. 
As a result, the cell is propelled at swimming speeds of 
about 30 /itm/s [1,3,5]. The process of bundling of nearby 
rotating "filaments" was studied in detail in macroscopic- 
scale experiments [7-9]. The cell tumbles and changes its 
direction of swimming randomly when one or more of the 
flagellar motors reverses its direction which forces the flag- 
ellar filaments to leave the bundle. In addition, a sequence 
of changes in the filament's handedness and pitch occurs 
[3-5,10]. Hence, the overall movement of a bacterium is 
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the result of altering intervals of tumbling and straight 
swimming. Chemotaxis steers the bacterium by just reg- 
ulating the tumbling frequency so that the net motion 
heads for a more favorable food environment [4,5]. 

For an object with a characteristic linear dimension 
a moving with velocity v through a Newtonian fluid, the 
ratio of inertial to viscous forces is given by the Reynolds 
number Re = avp/rj, where p is the fluid density and rj 
the viscosity [11,12]. Therefore, at low Reynolds numbers 
(Re <C 1), inertia does not play an important role, and the 
thrust pushing the object forward results solely from vis- 
cous drag. Swimming microorganisms in water are moving 
at very low Reynolds numbers [11]. E. coli bacteria, e.g., 
have a cell body of size a « 1/zm and move with veloci- 
ties of the order of v « 10/xm/s, which yields Re m 10~ 5 . 
Thus, the locomotion of microorganisms is fundamentally 
different from propulsion mechanisms in the macroscopic 
world (for comparison, a dolphin moves at Re rs 10 7 ). 

At low Reynolds numbers, the relative motion of two 
objects is governed by long-range hydrodynamic interac- 
tions which, to leading order, fall off with their inverse dis- 
tance [12]. They are also important in biological systems. 
Having in mind the propulsion mechanism of spermatozoa, 
Taylor modeled the hydrodynamics of two neighboring un- 
dulating tails, and found that hydrodynamic interactions 
synchronize the phases of lateral waves traveling down the 
tails [13]. Furthermore, the coordinated motion or stroke 
of beating cilia (known as metachronism) is believed to 
be mediated by hydrodynamic coupling [14-16]. In both 
cases, it is observed that the overall friction in the sys- 
tem is reduced by synchronization [13,16]. In analogy to 
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Fig. 2. The helix axis a, is 
tilted by the angle 0; against 
the z direction. The phase <j>i 
measures the angle between 
the rry-plane projection of the 
phase vector /3» and the x axis. 



to the anchoring strength of the terminal beads. Finally, 
we conclude in Sec. 5 discussing the role of the harmonic 
traps and the flexibility which they create. 



Fig. 1. Visualization of the helix geometry used in the sim- 
ulations (here with a phase difference of n/2). (a) All beads 
of one helix are connected ridigly with each other, (b) For the 
sake of clarity, the beads are "smeared" out along the helix. 
The top and bottom beads are anchored in harmonic traps. 
The illustrated helices are in their equilibrium positions (i.e., 
in the absence of driving torques). 

these examples, it was suggested that hydrodynamic in- 
teractions may also play an important role in how flagellar 
filaments synchronize their rotational motion so that they 
can form bundles [1] . 

In a recent paper, Kim and Powers studied hydrody- 
namic interactions between two rotating helices within the 
framework of slender-body theory [17]. The helices were 
considered as rigid and prevented from translation by ex- 
ternal forces, so that their axes were always parallel. The 
key result of this work was that there is no phase synchro- 
nization in this setup, when the two helices are driven with 
the same torque. 

In this paper, we consider a model which also consists 
of two stiff helices, thus neglecting any effects of elas- 
tic deformations. The helices are modeled by single beads 
that are rigidly connected to each other and are driven 
by constant and equal torques. In contrast to Rcf. [17], 
we "fix" the helices in space by anchoring their terminal 
beads in harmonic traps. This allows for slight shifts and 
tilts of the helices and thus implies some kind of flexibility, 
which is the major difference to Ref. [17]. In the following, 
we show that the phases of the rotating helices do indeed 
synchronize in this setup, and that the state of zero phase 
difference possesses lowest friction. 

The model is introduced in detail in Sec. 2. Then, sym- 
metry properties of the dynamic quantities are derived in 
Sec. 3. The numerical simulations of the helix dynamics 
are presented in detail in Sec. 4, where we analyze the 
data for phase synchronization, in particular with respect 



2 Model 

We consider two identical helices built of equal-sized beads 
[Fig. 1(a)] that are connected with each other by (virtual) 
rigid bonds. Thus, the helices cannot deform elastically. 
The centers of the beads are aligned along the backbone of 
the helix, with equal distances between successive beads. 

To describe the dynamics of the helices, we introduce 
body-fixed coordinate axes, given by the orthonormal vec- 
tors oti, (3i, and ai x /3j (i = 1, 2). The axis of a helix is 
represented by c*i, and the orientation of the perpendic- 
ular vector f3i shall describe the phase of the helix, i.e., 
the rotation about its own axis [Fig. 1(b)]. We define the 
phase angles 4n by the projection of (3i into the xy plane 
(Fig. 2). The angle between on and the z axis is the tilt 
angle 0j. 

The centers of mass of the helices are denoted by X{. 
The positions of the individual beads are then given by 

x v i = Xl + fte* + £ft + gat x 0i (1) 

with the internal coordinates 

h ( nm-l \ 

(2) 

to = r cos — v , 5, = r sin — v , 
m m 

where r is the radius of the helix and h its pitch. The 
bead index v runs from to nm — 1 for each helix, with m 
being the number of beads per winding and n the number 
of windings. 

The helices are driven by constant and equal torques 
that are always parallel to the respective helix axis, i.e., 
the torques are given by Doti with a fixed parameter D. 
Note that the assumption of a constant torque agrees with 
experimental studies of real flagellar motors [4,6], as we 
have already mentioned in our introductory remarks. To 
"fix" the helices in space, we attach single beads at the 
top and bottom end of each helix axis (Fig. 1) and anchor 
them in harmonic traps with equal force constants K. In 
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Fig. 3. Rotation of the two-helix 
system by 180° about the 2 axis (a) 
and about the x axis (b). The circu- 
lar arrows (indexed with i = 1, 2) on 
top of the helices denote that the re- 
spective helix is driven with torque 
Ti and rotates with velocity u>i in 
the indicated direction. The tubes 
are drawn as guide to the eye. 



equilibrium, both helix axes arc parallel, and their center- 
to-center distance is S. If one of the anchoring beads is 
displaced by Axf (where the index a refers to "top" or 
"bottom") relative to the center of the respective har- 
monic trap, the restoring single-particle force is 



F" = -KAx° 



(3) 



Finally, the total center- of-mass forces and torques acting 
on the rigid helices are 



Fi = Fl (with a = top, bottom) , 
Ti = Dot, + J2(xf - Xi) x Fl . 



(4) 



In the regime of low Reynolds numbers, the flow of 
an incompressible fluid with viscosity r\ obeys the quasi- 
static Stokes or creeping flow equations r/V 2 u — Vp = 
and V • u = [12, 18], where u is the flow field and p the 
hydrodynamic pressure. We assume the flow to vanish at 
infinity and impose stick boundary conditions on the sur- 
faces of all particles suspended in the fluid. The resulting 
flow field then couples the motion of the particles to each 
other. Due to the linearity of the Stokes equations, their 
translational and rotational velocities, Vi and depend 
linearly on all external forces and torques, Fj and Tj [12, 
18]: 

3 3 (5) 



is a 3 x 3 



Each of the mobilities nfj , fifj , nf. , and fi\ 3 
tensor, which couples the translations (superscript t) and 
rotations (superscript r) of particles i and j. They depend 
on the current spatial configuration of all suspended parti- 
cles. Since this dependence is highly nonlinear, they have 
to be calculated numerically. 

In our simulations, we use the numerical library hy- 
drolib [19] which yields the full set of mobility tensors 
for a given configuration of equal-sized spherical particles 
(based on the multipole expansion method). It implicitly 
accounts for (virtual) rigid bonds that keep the relative 
positions of the single beads in a rigid cluster fixed. Thus, 
hydrolib calculates an effective mobility matrix for the 
coupled center-of-mass translations and rotations, i.e., the 
indices i and j in Eq. (5) now refer to rigid clusters instead 
of individual beads (for details, see Ref. [19]). 



Therefore, with the forces and torques given in Eq. (4), 
we directly obtain the linear and angular velocities of the 
helices. The translational motion of the centers of mass is 
then governed by 

±i = Vi, (6) 

where the dot means time derivative. The rotational mo- 
tion of the helix axes oti and the phase vectors /3j follows 
from 

OLi = W, X OLi , 

(7) 

/3j = uji x fa . 

We integrate these equations in time by applying a second- 
order Rungc-Kutta scheme (also known as Heun algo- 
rithm) [20]. Note that the mobility matrices have to be 
evaluated at each time step since the positions and orien- 
tations of the helices change. 

While the trap constant K was varied to study the in- 
fluence of the anchoring strength on the helix dynamics, 
the driving torque D was kept fixed since it merely sets 
the time scale (given by the rotational frequency wo of an 
isolated helix). The time steps of the numerical integra- 
tion where chosen to correspond to roughly l/360th of a 
revolution of a single helix. 

The geometry of the two helices is shown in Fig. 1. 
Their backbones have a radius of r = 2.0a and a pitch 
of h = 6.0a, where a is the bead radius. The number of 
windings is n = 3, and the number of beads per winding 
is m — 5. The distance between the anchoring beads and 
the helix is the same as the pitch h. The equilibrium sep- 
aration of the helices, i.e., the distance of the upper/lower 
anchoring traps, is S — 7.0a. Note that the calculation of 
the mobility matrix is the most time consuming part in 
the simulations. Therefore, we had to restrict the number 
of beads in one helix. Furthermore, we will only present 
results for the set of parameters just introduced and con- 
centrate on the essential variable, namely the trap stiffness 
K. 



3 Symmetry considerations 

Consider for the moment two helices whose axes are com- 
pletely fixed in space, i.e., translation and tilt are pre- 
vented by appropriate forces and torques. In this case, 
the only remaining degrees of freedom are rotations about 
the axes of the helices. They are described by the phase 
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angles fa and the angular velocities fa — w*. We intro- 
duce the phase difference \ — fa — fa as the phase of 
the right helix relative to the left helix, when viewed as in 
Fig. 3(a,b) (left part). According to Eq. (5), the rotational 
velocities uui are functions of the phase angles fa since the 
mobility tensors depend on the spatial configuration. As 
the helices are driven by the same torques Tj = D about 
their axes, the synchronization rate is given by 

X = 0J 2 (fa,fa) - U 1 (fa,fa) = fJ,(fa,<fo)D, (8) 

where the effective mobility yu is 27r-periodic in fa. We 
choose this careful definition because we now want to de- 
rive symmetry properties of ji. 

We use the fact that the dynamics of the two-helix 
system must not change under arbitrary rotations of the 
whole geometry since the surrounding fluid is isotropic. By 
applying the two operations illustrated in Fig. 3, we create 
new configurations with left and right helices whose known 
dynamics we use to infer properties of the mobility /i. 

In the first case, we rotate the two-helix system by 180° 
about the z axis, as illustrated in Fig. 3(a). The velocities 
of the left and right helix are exchanged, i.e., toi <-> LU2 
and thus x — > ~X- On the other hand, the phase angles of 
the new left and right helix are, respectively, fa + tt and 
fa + 7T. Combining both statements, Eq. (8) yields 

fj,(fa +n,fa +tt) = -fj,(fa,fa) . (9) 

If the phases of the helices differ by tt (fa = fa + tt), one 
therefore obtains fi(fa,fa + tt) = —fi(fa,fa + tt) = or 

X = for x = tt, (10) 

i.e., the synchronization speed vanishes for any fa when- 
ever fa — fa — TT. 

Let us now rotate the two-helix system by 180° about 
the x axis [Fig. 3(b)]. Then the velocities of the left and 
right helix arc exchanged and reversed, i.e., lo\ <-> — oj 2 , 
and the synchronization speed x = — <*>i — (— ^2) = ^2 — wi 
stays the same. On the other hand, the angles transform 
as fa <-> — fa, and the driving torques are reversed, i.e., 
D — > — D. Again, combining both statements, Eq. (8) 
yields n(—fa, —fa)(—D) = fJ,(fa, fa)D, and thus 

-<f>i) = -(J-(<t>u<h)- (11) 

For helices of infinite length, the dynamics can only de- 
pend on the phase difference x and not on the single 
phases fa . This is obvious since a phase shift of both he- 
lices is equivalent to a translation along the helix axes, 
which does not change the dynamics. Hence, Eq. (11) 
reads /u(x) = — m(x) = 0> i- c -i f° r parallel helices of in- 
finite length, the synchronization rate x = ^{x)D van- 
ishes for any phase difference x, and therefore, they do 
not synchronize towards x = [21]. 

4 Synchronization 

We now study the rotational dynamics of two helices whose 
terminal beads are anchored in harmonic traps of finite 
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Fig. 4. Synchronization of the helix rotations. The phase dif- 
ference of the two helices tends towards zero, starting from x 
slightly smaller than tt. The symbols are simulation data at 
two different trap strengths (values of K in units of D/a). For 
clarity, not every data point is plotted. The solid line shows 
the master curve of Eq. (12). The insets enlarge the small os- 
cillations at x ~ and tt (here for the case K — 0.05 D/a, but 
the amplitudes do not depend strongly on K). Note that the 
scaling of the two insets is the same. 

strength, as introduced in Sec. 2. Thus, the helices can 
be shifted and tilted, and their axes undergo a precession- 
like motion while each helix itself is rotating about its 
respective axis. The orientations of the helices in space 
are described by the vectors on and /3j [see Fig. 1(b)] and 
the corresponding angles 9i and fa, as defined in Fig. 2. 

Figure 4 shows the phase difference x^^ — fa for two 
trap stiffnesses K as a function of a reduced time t(K), to 
be defined below. Starting with x slightly smaller than tt, 
the phase difference decreases continuously (with steep- 
est slope at x = 7r /2) and finally approaches zero, i.e., 
the two helices do indeed synchronize their phases. The 
simulations reveal that the dynamics does not depend sig- 
nificantly on the phases fa themselves, but is predomi- 
nantly determined by the phase difference x- Note that 
this feature may be expected since the dynamics of par- 
allel helices of infinite length can only depend on x as 
explained in Sec. 3. To be concrete, we observe that the 
rotational velocities fa undergo oscillations of only about 
1% around a mean value during one rotational period. 
Their amplitude does not depend on the trap stiffness K, 
i.e., the oscillations originate from the slight dependence 
on the phases themselves (and not from the precession of 
the axes). Therefore, starting with different values for fa 
but the same value of x yields the same curve (except for 
differences in the small oscillations illustrated in the insets 
of Fig. 4). 

The mean rotational velocities, averaged over one rota- 
tional period, increase during the synchronization process 
from about 0.92ojo at x ~ ^ to about 0.95wo at x ~ 
(Fig. 5). Thus, the hydrodynamic drag acting on the he- 
lices is minimized during phase synchronization (see also 
Ref. [17]). Since the torques are constant, the dissipation 
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Fig. 5. Mean rotational velocities, averaged over one revolu- 
tion (in units of the rotational frequency loq of an isolated he- 
lix). The oscillations about the mean value are of the order of a 
few percent, they decrease slightly during the synchronization 
process. The example shown is for trap strength K = 0.15 D/a. 
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Fig. 6. Synchronization speed (taken at a relative phase \ w 
7r/2, as illustrated in the inset) as a function of inverse trap 
strength K~ x . The frequency scale ojo for the synchronization 
speed is the angular velocity of an isolated helix. The symbols 
indicate values extracted from simulations at different K. The 
dashed line is an empirical fit (see text) . 



rate Tj-Wj is maximized. This observations agrees with 
the interesting fact that the Stokes equations can be de- 
rived from a variational principle where one searches for 
an extremum of the dissipated energy J <JijAijd 3 r (mj is 
the stress tensor and Aij the symmetrized velocity gradi- 
ent) under the constraint that the fluid is incompressible 
[22,23]. The pressure enters via the Lagrange parameter 
associated with the constraint. 

In Sec. 3, we showed for fixed parallel helices, based 
on pure symmetry arguments, that their synchronization 
rate vanishes for a phase difference of x — tt [see Eq. (10)]. 
In that case, the two-helix configuration is symmetric with 
respect to a rotation by 180° about the z axis [see Fig. 3(a)] . 
Our reasoning of Sec. 3 can be extended to the case of 
non-parallel helix clXGS ; clS long as the same symmetry is 
preserved. However, x — n does not correspond to a stable 
state. Starting with \ marginally smaller than it, the sys- 
tem tends towards phase difference zero. The simulation 
with K = 0.05 in Fig. 4 was launched, e.g., at \ — 0.9947T 
with both helices in equilibrium position and orientation, 
as shown in Fig. 1. 

On the other hand, the synchronized state \ — is sta- 
ble against small perturbations since configurations with 
X between and — 7r synchronize towards zero phase differ- 
ence, too. This was checked by simulations, but can also 
be derived from Eq. (9). The corresponding rotation of 
Fig. 3(a) creates new left and right helices with a change 
in sign for \ an d X relative to the original helices which 
explains our statement. Furthermore, starting a simula- 
tion with exactly x = 0, the helices remain synchronized 
on average (i.e., (x) = 0), but there are still small oscil- 
lations as illustrated in the lower right inset in Fig. 4 for 
the case where both helices started in equilibrium position 
and orientation. 

Averaging over small oscillations, we find that the re- 
sulting smoothed curves for the phase difference \ obey 



an empirical law of the form 



x(t) = -(1-tanhr), 



where 



r(K) = -(t- t n/2 ) 



At 



(12) 



(13) 



is the reduced time, already mentioned above, and t n de- 
notes the time where x = 7T /2, i.e., the location of the in- 
flection point. Its slope \dx/dt\t=t„ /2 depends on the trap 
stiffness K and so does r. As Fig. 4 strikingly reveals, this 
law works very well. By plotting the phase difference x 
versus the reduced time t(K), the curves collapse on the 
master curve given by Eq. (12). Since the dynamics at low 
Reynolds numbers is completely overdamped, we expect 
this law to follow from a differential equation which is of 
first order in time. Taking the first derivative of Eq. (12) 
with respect to r, we find that x( T ) obeys the nonlin- 
ear equation x(t) = (2/tt)x(t)[tt — x( T )L known as the 
Verhulst equation and originally proposed to model the 
development of a breeding population [24] . However, it is 
not clear how to derive this equation from first principles 
in our case. 

An important result is that the speed of the synchro- 
nization process decreases with increasing trap stiffness 
K . The values plotted in Fig. 6 for different K~ x are the 
slopes |dx/di|t=t 7r/2 extracted from simulation data at the 
inflection point with a relative phase of x ~ 7r /2- The 
curve in Fig. 6 can be extrapolated by the analytic form 
ci tanh c 2 K~ 1 (dashed curve), where the fit parameters as- 
sume the values c\ = 3.67 • 10~ 3 wo and c 2 = 0.0685 D/a. 
In the limit of infinite trap strength, i.e., for K^ 1 — > 0, 
the synchronization speed clearly tends towards zero, i.e., 
an infinitely strong anchoring of the helix axes does not 
allow for phase synchronization. 
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Fig. 7. Tilt of the left and right helix axis for K = 0.05 and 
0.15 as a function of reduced time r. The fast periodic oscil- 
lations of the tilt angles 9i yield the black bands since they 
cannot be resolved on the time scale used here. The insets vi- 
sualize the precession of the helix axes by showing the tip of 
the vectors q, (at \ ~ t/2)- From the outer to the inner "cir- 
cular" orbit, the trap stiffness assumes the values K — 0.05, 
0.1, 0.15, and 1.0 (in units of D/a). 



In Fig. 7, we illustrate how the tilt angles Qi (for their 
definition, see Fig. 2) vary during the synchronization pro- 
cess. The mean tilt angle as well as the amplitude of its 
periodic oscillations decrease when the phase difference 
approaches zero. Obviously, the dynamics of the helices 
depends on the stiffness of the harmonic anchoring of the 
top and bottom terminal beads. In a weaker trap, the tilt 
of the helix axes out of equilibrium is more pronounced 
compared to a stronger trap. The insets in Fig. 7 track 
the precession-like motions of the helix axes. The stronger 
the trap, the smaller the radius of the "orbit" or the tilt 
angle. [Since the simulations were started with both axes 
aligned along their equilibrium direction, the trajectories 
(oti X , a iy ) first move radially away from the origin and then 
enter the "precession orbit" .] 

In real flagellar motors, the torques on the two heli- 
cal filaments are not exactly the same. To test whether 
the phenomenon of synchronization still occurs, we now 
consider slightly different driving torques for the two he- 
lices. The first helix is still driven with torque D, while the 
second helix is driven with D + AD (AD > 0). We ob- 
serve that, for torque differences AD below a critical value 
AD C , the two helices indeed synchronize towards a phase 
difference Xoo = x(i — > oo), which, in general, is not zero. 
The results are shown in Fig. 8, where we plot 1 — Xoo/f 
as a function of AD/AD C . With increasing AD, the phase 
lag Xoo increases from zero to tt/2. For AD > AD C , there 
is no synchronization, and the phase difference grows con- 
tinuously. Note that the reduced critical torque difference 
AD C /D = 3.22- 1CT 3 corresponds to the reduced frequency 
difference x/luo — (lo2 — u>i)/u>o observed at x — 7r /2 for 
equal torques (AD = 0). At the critical torque difference, 
the helices synchronize towards a phase lag of \ = tt/2. 
This means that AD C just compensates the difference in 



i 




0.6 0.8 

AD/AD C 

Fig. 8. Phase difference Xoo (plotted as 1 — Xoo/f ) of the syn- 
chronized state as a function of the torque difference AD. The 
phase difference increases from Xoo = to tt/2 (or 1 — Xoo/f 
decreases from 1 to 0) with increasing torque difference. For 
AD > AD C , there is no synchronization. 



the effective mobilities of the two helices, which is largest 
for x = 7r /2; thus the helices rotate with the same speed. 

At the end, we mention that all results presented here 
refer to helices whose rotational direction is given in Fig. 
1(b). Reversing the direction of rotation does not change 
the dynamics of the two-helix system since this can also 
be achieved by the operation shown in Fig. 3(b) that does 
not change the synchronization speed. 



5 Conclusions 

We have reported that two rigid helices whose terminal 
beads are anchored in harmonic traps and which are driven 
by equal torques synchronize to zero phase difference. The 
effect is robust, i.e., if the torques are unequal, the helices 
synchronize to a non-zero phase lag below a critical torque 
difference. This agrees with observations in Ref. [8] where 
the helices do not bundle if the motor speeds are suffi- 
ciently different. Increasing the stiffness of the anchoring 
traps, decreases the synchronization rate. We attribute 
this to the jiggling motion of the two helix axes which is 
more and more restrained. 

In the limit of infinite trap strength, our results are 
consistent with recent work based on slender-body the- 
ory for two rigid helices [17]. If the helices are prevented 
from translation and their axes are always kept parallel, 
then there is no synchronization possible. Therefore, we 
conclude that the additional degree of freedom due to the 
finite anchoring of the helix axis, i.e., the jiggling motion, 
is essential to enable phase synchronization in our model. 

At a first glance, our model might appear too artifi- 
cial for describing the hydrodynamic coupling of flagella. 
However, our results clearly indicate that some kind of 
flexibility is essential to allow for phase synchronization. 
In reality, this flexibility might have its origin in elastic 
deformations of the rotating flagella. Therefore, the next 
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step is to make the flexural and torsional stiffness of the 
helices in our model finite. 

The helix used in our numerical investigation with ra- 
dius r — 2a and pitch h = 6a is far from the values of a 
real bacterial filament with r = 20a and h = 200a where 
a is now half the filament diameter. We have tried differ- 
ent paths in the parameter space (r, h) to connect both 
cases, i.e, our "fat" helix and the real slender helix. We 
observe that making the helix more slender decreases the 
synchronization speed in units of u>o- This makes sense 
since the induced flow from the rotation of slender he- 
lices is smaller. To reduce computer time, we extrapolated 
the synchronization speeds for different paths in the (r, h) 
space towards the real helical flagellum and found that the 
speed is reduced by a factor of 60 to 70 compared to the 
results reported in this article. At a first glance this seems 
to be a discouraging result. However, from preliminary re- 
sults of helices with finite bending and torsional flexibility, 
we know already that these factors considerably increase 
the synchronization speed so that it will be of biological 
relevance. 

Furthermore, we made a comparison of the resistance 
matrix of a real flagellum modeled by a sequence of spheres 
with resistive force theory as summarized, e.g., in Rcf. [25] . 
We found that for motion and rotation along the helix 
axis, the single matrix elements differ by less than a factor 
of two. This convinces us that modelling a flagellum with 
the method presented in this article is appropriate. For 
our fat helix, however, the deviations are larger since the 
filament is relatively thick compared to radius and pitch 
of the helix. So the conditions for the validity of resistive 
force theory are not satisfied so well. 

We also checked whether it is important if the helices 
are forced to stay at their position or if they are allowed 
to propel themselves. This was done by letting the helices 
move along the z axis but still keeping them in harmonic 
traps along the x and y direction. However, we did not 
see a significant difference in the dynamics of the synchro- 
nization process compared to the results presented in this 
work. 

As a final remark, we point out that the following gen- 
eral mechanism may exist in systems of low Reynolds num- 
bers: Synchronizing the motions of some objects via hy- 
drodynamic interactions needs some kind of "flexibility" . 
If the motions of the objects are constrained too much, 
synchronization cannot occur. We have observed a simi- 
lar behavior for particles circling in a toroidal harmonic 
trap and driven by a constant tangential force [26]. After 
some transition regime, the particles reach a synchronized 
state where they perform a periodic limit cycle. We ob- 
serve that for increasing trap stiffness, i.e., for decreasing 
oscillations along the radial direction, the time to reach 
this limit cycle increases. 

This work was supported by the Deutsche Forschungsgemein- 
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